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ABSTRACT 

When an image of a strongly lensed quasar is microlensed, the different components of its spectrum are expected to be differentially 
magnified owing to the different sizes of the corresponding emitting region. Chromatic changes are expected to be observed in the 
continuum while the emission lines should be deformed as a function of the size, geometry and kinematics of the regions from which 
they originate. Microlensing of the emission lines has been reported only in a handful of systems so far. In this paper we search 
for microlensing deformations of the optical spectra of pairs of images in 17 lensed quasars with bolometric luminosities between 
2Q44.7-47.4 er g/- s an( j black hole masses 10 76 ~ 9 8 Mq. This sample is composed of 13 pairs of previously unpublished spectra and four 
pairs of spectra from literature. Our analysis is based on a simple spectral decomposition technique which allows us to isolate the 
microlensed fraction of the flux independently of a detailed modeling of the quasar emission lines. Using this technique, we detect 
microlensing of the continuum in 85% of the systems. Among them, 80% show microlensing of the broad emission lines. Focusing 
on the most common emission lines in our spectra (C III] and Mg II) we detect microlensing of either the blue or the red wing, or 
of both wings with the same amplitude. This observation implies that the broad line region is not in general spherically symmetric. 
In addition, the frequent detection of microlensing of the blue and red wings independently but not simultaneously with a different 
amplitude, does not support existing microlensing simulations of a biconical outflow. Our analysis also provides the intrinsic flux 
ratio between the lensed images and the magnitude of the microlensing affecting the continuum. These two quantities are particularly 
relevant for the determination of the fraction of matter in clumpy form in galaxies and for the detection of dark matter substructures 
via the identification of flux ratio anomalies. 
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1. Introduction 

Soon after the discovery of the first gravitationally lensed quasar, 
it has been realised that microlensing (ML) produced by com- 
pact objects in the lensing galaxy towards a multiply imaged 
quasar could be used as an astrophysical tool to probe the in- 
ner parsecs of distant quasar s (|Chang &T<e fsdal 1979{ [Kayser| 
[eTaLl[T98r3l |Paczynski|[T98r5} |Grieger et al.||1988) . Microlenses 
typically magnify regions of the source on scales similar to 
or smaller than a few micro-arcsecs, the size of their angu- 



lar Einstein radius R E (Wambsganss 1998 2006 Schmidt & 
Wambsganss 2010). Hence, the quasar continuum region (accre- 



tion disc) and the broad line region are likely to be microlensed. 
Because the magnification varies with the source size, the power 
law continuum emission of quasars is expected to be more mag- 
nified as the wavelength, and hence the continuum size, de- 
creases. Microlensing is therefore expected to produce signifi- 
cant color changes in macro-lensed quasar images (Wambsganss 



* Based on observations made with the ESO-VLT Unit Telescope 
# 2 Kueyen (Cerro Paranal, Chile; Proposals 074.A-0563, 075.A-0377, 



077.A-0155, PI: G . Meylan). 



** The new spectra presented in this paper are available in elec- 
tronic form at the CDS via anonymous ftp to cdsarc.u-strasbg.fr| 



(130.79.128.5) or via http : //cdsweb . u- strasbg . fr/cgi-bin/ 



qcat?J/A+A/ and via the German virtual observatory http://dc. 
g- vo . org/mlqso/q/web/form 



& Paczyriski|199i). T he latter have indeed been observed (e.g. 
Wisotzki et al.|19931|CTa"eskens et al.|[200T] |Burud etlLl[2002l 
Wucknitz et al.|2003] ). The effect on the broad line emitting re- 
gion (BLR) has been first addressed by Nemiroff (1988) who 
calculated the changes produced by a single microlensing-star 
on the emission lines. This analysis has been refined soon after 
by Schneider & Wambsganss ( 1990| ) who considered the more 
realistic case of ML by a population of microlenses. These pa- 
pers have demonstrated that microlensing of the BLR could be 
significant and does depend on the BLR geometry and its kine- 
matics. They also showed that ML of a spherically symmetric 
BLR (in geometry and velocity field) would lead to symmetric 
variations of the emission lines (i.e. of the blue and red com- 
ponents) while ML of a keplerian disc would lead in general to 
asymmetric variations of the emission lines and possible shift 
of the line centroid. Microlensing affecting more peculiar line 
profiles from e.g. broad absorption lines quasars, or generated 
in a relativistic disc have been discussed i n |Hutsemekers | (| 1 99 3 \ ; 
Hutsemek ers et al.|(fl994l);|Lewis & Belle| ( |1998| ); |Belle & Lewis| 
(2000); Popov icet al.| ( |2001[ ). Despite these promising results, 
detection of microlensing in the emission lines remained elusive 
(Filippenko 1989 Lewis et al. 1998 1 or invoked to explain differ- 



ences between the spectra of candidate lensed quasars ( Steidel & 
|Sargent|199~T)|Small et al.|1997| ). T he interest in BLR mic rolens 



ing got revived with the papers of lAbajas et al. ( 2002 ) and of 



Lewis & Ibata (2004), who re-investigated this question after 
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the discovery that the BLR was smaller than previously thought 
( |Wandel et aL][T999l |Kaspi et al.||2000] >. Based on the size of 
the BLR me asured in NG C5548 and using the scaling relation 
Rblr °c L 7 , |Abajas et al.| ( |2002| estimated that the BLR should 
be significantly microlensed in about ~ 30% of the systems. 
They also extended the work of Nemiroff ( 1988), and calculated 
the microlensing by a single lens for various BLR geometries 



and kinematics, considering BLR models described in Robinson 



(1995 



(2002 



micro 



et al. 



Lewis & Ibata (2004) extended the work of Abajas et al. 



by using more realistic microlensing patterns. Finally, 



ensing of a biconical BLR, already presented in Abajas 
(2002 1 and in Lewis & Ibata (2004), for two peculiar ori 



entations of the axis of the bicone, has been discussed for more 



general bicone configurations in Abajas et al. (2007 1. These pa- 
pers have confirmed most of the earlier findings and made more 
detailed predictions on the line-shifts and asymetries induced by 
ML. They also showed that the line deformation depends only 
weakly on the value of the surface density k and shear y at the 
position of the lensed images, but more strongly on the orienta- 
tion on the BLR w.r.t. to the direction of the shear. 

Many papers dedicated to the detection and interpretation 
of a microlensing signal have focused on the Einstein Cross 
= Q2237+0305, which is probably the most favourable object 
for microlensing studies. Indeed, this system has a negligible 
time-delay between the lensed images, which enables one to 
easily disentangle microlensing and intrinsic variability, and a 
low-redshift lensing galaxy which leads to a small Re and to 
relatively fast microlensing variations (Mosquera & Kochanek 
201 1) . After the first detection of microlensing in the contin- 



uum and in the broad line ( |Irwin et al.|1989| |Filippenko|1 989), 
microlensing has started to be used as a tool to constrain the 
size of the accretion disc and of the broad fine emitting region 



jLewis et al |1998[|Wyifhe et al.|2000b"l|Yonehara|200T]|Wyithe 
|et al.|2000a) . In the last decade, important progress in observa- 
tional techniques allowed photometric and spectrophotometric 
monitoring of the individual lensed images to be successfully 
carried out (Wozniak et al.|20 00; Dai et al. 2003; Anguita et al. 
2008b] |Eigenbrod et al.|2 008 ; Zim meret al'|2011| >. On the other 



hand, the development of more advanced numerical techniques 
allowed quicker calculation of magnification maps and more so- 
phisticated analysis (Kochanek 2004 Poindexter & Kochanek 
2010 a~lb"l|B"ate et al.|2010||Mediavilla et al.|2011t|Garsden et al. 
201 1 ||Bate & FlukepoT^ . Owing to these two ingredients, tight 



constraints on the size of the accretion disc and on its tempera 
ture profile, on the size of the broad line region and on its kine- 
matics have been obtained for Q22 37+0305 fEigenbrod et al. 
2008j|Poindexter & Kochanek|2010a||O'Dowd et al |201 lHSluse 
etal.|2011) . ~ 



Most of the recent papers focused on the study of the quasar 
accretion disc, which can be done using broad band photometry 
from X-ray to optical and near-infrared wavelengths (e.gJPooley 
et al.|2007||Floyd et al.|2009l|Hutsemekers et al.|20To||Dai et aT 
2010] |Blackburne et al.||2011| |Mufioz et al.||2011| > Studies of 
BLR microlensing are more sparse and detections have been re- 
ported only for a handful of systems ( Richards et al.|2004[|Wayth 



et al.|2005HKeeton et al.||20D6| |Sluse et al.|2007HHutsemekers 



et al.|2010| . In this paper, we present a careful re-extraction and 
analysis of archive spectra of a sample of 13 lensed quasars ini- 
tially observed with the aim of measuring the redshift of the lens- 
ing galaxy. Our systematic analysis allow us to characterise the 
microlensing-induced deformation of the emission lines. To get 
a more complete overview of the microlensing signal, we also 
discuss the signal detected in four objects we presented else- 
where. 



From our spectra, we also derive flux ratios corrected for mi- 
crolensing which are closer to the intrinsic flux ratios between 
pairs of images. This is important for the study of doubly im- 
aged quasars for which the flux ratios are mandatory to con- 
strain the lens models, because of the few observational con- 
straints in these systems (e.g. Chantry et al.||20i0 Sluse et al. 
|2012[ ). Intrinsic flux ratios are also particularly relevant for the 
identification of flux ratio anomalies possibly produced by dark 
matter substructures or ML. A popular technique to study the 
amount of massive dark matter substructures in galaxies relies 
on the identification of flux ratios between lensed image pairs 
which deviate from lens model prediction, i.e. the so called flux 
ratio anomaly (|Mao & Schneider|1998| | Metcarf & Madau|2001[ 



Dalai & Kochanek||2002| |Keeton et al.|2003| |Fadely & Keeton 
|201 1; Zackriss on & Riehm|2010| and references therein). One 
of the current hmitation of this technique is the small number of 
reliable flux ratios which may be used to identify an anomaly. 
Indeed, most of the lensed quasars are observed at visible and 
near-infrared wavelengths where microlensing and differential 
extinction significantly contaminate the flux ratios, while only a 
handful of systems are detected in the mid-infrared or at radio 
wavelengths where these two effects are negligible. In this paper 
we discuss how flux ratios derived using spectra may provide a 
good proxy to intrinsic flux ratios, allowing one to significantly 
extend the sample of objects where flux ratio anomalies can be 
studied. 

The structure of our paper is as follows. In Sect.|2j we present 
the extraction and flux calibration of the archive spectra, and the 
spectra from literature. In Sect. [3] we derive the physical prop- 
erties (bolometric luminosity, black hole mass, Eddington ratio) 
of the lensed quasars. We also discuss there the published red- 
shifts of the quasar and provide an alternative value based on 
the Mgll emission line. Sect.|4]is devoted to the analysis of the 
microlensing in the spectra. It includes a presentation of the tech- 
nique used to isolate the microlensed fraction of the flux, a de- 
scription of the microlensing signal observed in each object and 
a discussion on the accuracy of our microlensing-corrected flux 
ratios. In Sect. [5 we discuss the microlensing signal observed in 
the continuum, the occurrence and variety of microlensing of the 
broad lines and the consequences for the structure of the BLR. 
Finally, we summarize our main results in Sect. [6] 



2. Observations and data processing 

2.1. Observations 

We have gathered spectroscopic observations of 13 gravitation- 
ally lensed quasars previously presented in |Eigenbrod et al. 
(2006b 2007). The spectra were obtained in multi-object mode 
with the FORS1 instrument mounted on the ESO-VLT Unit 
Telescope 2 Kueyen (Cerro Paranal, Chile). The object names 
and main observational characteristics of the data are presented 
in Table [T] The following information is provided: name of the 
object, observing date, total exposure time, average seeing, aver- 
age airmass, and reference where more details regarding the ob- 
servations and data reduction can be found. All the spectra were 
obtained through 1" slit width with the G300V grism and the 
GG435 order sorting filter. This setup provides us with spectra 
from 4400 to 8650 A. All the objects, except SDSS J1226-0006, 
were observed with the High Resolution collimator of FORS, 
leading to a spectral resolution R = A/AA ~ 210 at 5900 A, and 
to a pixel scale of 2.9 A in the spectral direction and of 0. 1 " in the 
spatial direction. For SDSS J1226-006, the Standard Resolution 
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Notes: f New measurements, see Sect. |3.4| 

References. (l) |Eigenbrod"ltlil ( |2006b) ; (2) |Eigenbrod"eraLlp007l > 

Table 1. Data summary for our sample of 13 lenses. We give the object name, the redshift of the QSO lensed images (z s ), of the 
lensing galaxy (zi), the names of the lensed images in the slit, the date of observation and exposure time, average seeing, airmass 
and reference to the paper where the data where first presented. 



mode of the FORS instrument was used, leading to 0.2" per pixel 
in the spatial direction and to R ~ 400 at 5900 A. 

2.2. Reduction and deconvolution 

We used the combined 2D spectra of each object as obtained 
after the standard reduction procedure presented in Eigenbrod 
|et al.| ( |2006a| |2007[ l. An example of such a 2-D spectrum is 
displayed in FigTjTji. The extraction of the spectra and the flux 
calibration, as described hereafter, have been updated compared 
to the original publication. Note that the differences arise be- 
cause the scientific goal of the original papers was the determi- 
nation of the redshift of the lensing galaxy, while we are inter- 
ested here in a detailed study of the unpublished quasar spec- 
tra. Like |Eigenbrod et al.| ( |2006a| |2007|>, we use MCS dec on- 
volution algorithm adapted to spectra ( |Courbin et al. 2000]) in 
order to deblend the components of the lens system. This al- 
gorithm deconvolves the observed frame into a frame of finite 
resolution with a Gaussian PSF chosen to have FWHM = 2 pix- 
els. During the deconvolution process, the flux is separated in 2 
channels, one channel (i.e. "the point-source channel") contain- 
ing only the deconvolved flux of the point-like sources, and one 
numerical channel (i.e. the "extended channel", Fig.[T]:) contain- 
ing the remaining flux correlated over several pixels. This ex- 
tended channel contains mostly signal from the lensing galaxy 
and sometimes improperly subtracted sky signal. We improved 
the extraction of the spectra with respect to the former decon- 
volution by systematically assuming 3 point sources in the 2D 
deconvolved spectra (Fig. [TJ)): two for the QSO lensed images 
and one for the lens galaxy which has a relatively peaked pro- 
file in the center. Flux from the galaxy which deviates from a 
PSF is included in the "extended channel". The residuals are 
systematically inspected (Fig. [TJl) to identify possible deconvo- 
lution problems. For objects where the deconvolution was not 
satisfactory, we also tested alternative PSFs and chose the one 
leading to the lower residuals. Finally, the 1-D spectra of the 
QSO and of the lensing galaxy (Fig. [T£) are corrected from 



the response curve based on standard stars, and from differen- 
tial extinction using the updated Paranal extinction published by 
Patat et al. (201 1). An example of 1-D flux calibrated spectra for 
SDSS J1335+0118 is displayed in Fig. [J, g, h. The extracted 
spectra are available on electronic form via CDS, while 1-D and 
2-D spectra are available via the German Virtual Observatory^] 

When possible, we checked that the flux ratios derived from 
our spectra are compatible with nearly simultaneous R-band 
ratios from literature or from the COSMOGRAIL monitoring 
project. We could not obtain simultaneous R-band flux ratios for 
the following systems: SDSS J0246-0825, SDSS J0806+2006, 
FBQ 0951+2635, BRI 0952-0115, SDSS J1138+0314 and 
Q 1355-2257. Photometric and spectro-photometric data always 
agree within the error bars except for WFI 2033-4723 for which 
we find C/B ~ 0.85 while |Vuissoz et al] (f2008|, retrieved 
C/B = 0.7 in the r-band on about the same date (MJD=53500). 
Because of that disagreement, we have multiplied the spectrum 
of image B in WFI 2033-4723 by 1.2 in order to match the pho- 
tometric flux ratios. Our comparison with photometric flux ratios 
implies that the systematic error on the spectral ratio introduced 
by slit losses is < 5 %. 

2.3. The extended sample 

In order to provide a more global overview of the variety of 
signals which may be observed, we completed our main sam- 
ple with five systems we studied elsewhere: SDSS J0924+0219 
(Eigenbrod et al.|2006a|l, H14134-1 17 (|Hutsemekers et al.|2010 



J1131 -1231 ( |Sluse et al.||2007]>, HE 2149-2745 (|Burud et a 
20021), and Q2237+0305 ( |Eigenbrod et~aLl|2008| |Sluse et a 



|201 1[ ). This is not a complete list of known microlensed quasars 
(see e.g. |Richards et~aT1|2004[ |Anguita et iLl|2008a| ), but these 
systems were among the first quasars lensed by a single lens- 
ing galaxy where significant microlensing of the BLR has been 
unveiled and discussed, with the advantage that the fully pro- 
cessed spectra are at our disposal. In the following, we refer to 
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Fig. 1. Example of deconvolution for SDSS J1335+01 18. The first two lines show the original frame (a) and its deconvolved version 
spatially oversampled by a factor 2 (b), the deconvolved flux which is not PSF-like (c) and the residual frame (d) with upper/lower 
cuts at +3<x. The location of the centroid of the deconvolved PSFs is indicated with solid lines in panels (b, c, d). Panel (e) shows 
the extracted spectra of the lensed QSO images and of the lensing galaxy derived from the deconvolution. The last three panels (f, 
g, h) show the calibrated spectra of these components following the procedure explained in Sect. 



2.2 



this sample as the "extended sample". In order to derive the lumi- 
nosity of Jl 131-123 1 and H1413+117 as for the other systems, 
we performed an approximate flux calibration of their spectra, 
by matching the broad-band magnitudes derived from the spec- 
tra with published absolute photometry. We used the /-band 
magnitude of |Sluse et al^pOOrj) for Jl 131-1231 and the //-band 
2MASS photometry Z/=14.531+0.054 for H1413+117. 



3. Physical properties of the sample 

In this section, we calculate the black hole mass, the absolute 
luminosity, and the BLR size of each of the lensed quasars. We 
explain in Sect.|3.1||3.2| & 3.3 how we derive these quantities. In 



Sect. 3.4 we present updated redshifts of the quasars based on 
the Mg II line. Finally, we present the results of our calculation 
and the error estimates in Sect. 3.5 We emphasize that the de- 
termination of the luminosity makes use of macro-magnification 
ratios presented in Sect. [4] 



3.1. Black Hole mass measurement 

We estimate the virial black-hole mass using the width of 
the Mgll emission line, following the original prescription of 
|McLure & Dunlop ( |2004| l. In order to convert line width and 
continuum luminosity into black hole masses, we use a relation 
of the form 



log 



/ M BH \ 

\io 6 mJ 



a +/31og 



-3000 



10 44 erg/s 



-ylog 



/ FWHM ^ 
\ 1000km/sy 



(1) 



where Mbh is the black hole mass, L3000 = A/a(3000A) in 
erg/s and FWHM is the Full Width at Half Maximum of the 
Mgll emission line. We use a = 1.13 ± 0.27, /3 = 0.5 and 
y = 1.51 ± 0.49 as derived by |Wang et aT] ( [20091 E q- 10 ) from 
an analysis of a sample of about 500 SDSS spectra of intermedi- 
ate redshift quasars. We derive L3000 an d FWHM(Mg II) by si- 
multaneously fitting two Gaussians to the Mgll line on top of a 
pseudo-continuum component. The latter is the sum of a power- 
law continuum and of a broadened Fell template obtained by 
convolving (in the velocity space) the Fell template derived 
from I Zw 1 by Vestergaard & Wilkes (2001 ), with a Gaussian 
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of adequate width to reproduce the observed Fe II emission. The 
fit is always performed in the range [2200, 2675]A , excluding 
regions affected by the sky or intervening absorption. In order to 
use a methodology similar to the one commonly used, and de- 
spite the fact that the evidence for a narrow Mg II component 
in AGNs remains debated ( e.g. |McLure & Dunl op 2004; Sluse 
|et al.|2007||Fine et al |200"8"t|Sheti et al.|2008[|Wang et al.|2 009), 
we model the Mgll emission with a narrow (FWHM < 1200 
km/s) and a broad component {FWHM > 1200 km/s). Such a 
fitting procedure is similar to those used by McLure & Dunlop 



(2004} and|Shen et al.| (2 008 > but differs slightly from the one 



used by Wang et al. (2009) who included a Balmer continuum 
component in their fit and explicitly model Mg II as a doublet. 
Because of the low resolution of our spectra, it is not meaning- 
ful to split the Mg II /12798 A line into its individual components 
/12796//12803 which are spectroscopically unresolved. On the 
other hand, we have not fitted a Balmer continuum because of 
the degeneracy between this component and the power law in- 
dex. Not accounting for the Balmer continuum will not affect 
the measurement of the FWHM but could lead to an overesti- 
mate of ^3000- The Balmer continuum is commonly assumed to 
contribute to ~ 10% of L3000, but recent results suggest a contri- 
bution of ~25% in average, with a large scatter (Wang, private 
communication, 



enough to weak 



Jin et al.| 201 1 1. This remains in general small 



y affect Mbh since Mbh 00 ^3000- During the 
fit, we test for three different Fe II template^] (a) a "low level" 
template which assumes that there is no Fell emission under 
Mg II, (b) a "fiducial" template which assumes that the intensity 
of Fell in the range ~ [2757, 2825] A equals the ave rage inten- 
sity of Fell measured in the range [2930, 2970] A (Fine et al. 
2008), (c) a "high level" template which assumes that the inten- 
sity of Fell in the range ~ [2757, 2825] A is 40% larger than 
in case (b), such that the intensity of Fell under Mgll equals 
the amount of Fe II on each side of the line. We use the fidu- 
cial template to estimate FWHM(Mg II). Because the fraction 
of Fe II under Mg II is the main source of error in the measure- 
ment of the FWHM, we use the spread between (a), (b) and (c) 
to set (Tfwhm- The luminosity L3000 is estimated from the power 
law continuum and corrected from the macro-magnification as- 
sociated to lensing (see Sect. 3.2 1. 



In four objects, we were unable to measure the black hole 
mass from Mgll, because this line falls outside of the wave- 
length range covered by our spectra or is located less than 
300 A from the edge of the spectra (HE 2149-2745). For these 
systems we use black hole mass estimates from literature ob- 
tained using a relation similar to Eq{T] but for the C IV line. 
Black-hole mass estimates derived from that line are more 
controversial because of uncertainties on the geometry of the 
C IV emitting region, because of the contamination from a pu- 
tative narrow CIV emission and because virial equilibrium is 
not expected to hold perfectly for such high ionization lines 
~ 2002) |Bachev et al.||2004l |Marziani et~aL] 



( Rich ards et al 
2006 1 Fine et al. 



|20TT] |Assef et al.||20TT |Wang et al.||20TT 
Richards et al.||20lT Marziani & Sulentic 2012 1. For two sys 
terns (H1413+117 and Jl 131-1231), we also used H/3 to calcu 
late Mbh- For this purpose, we used the relation (A7) of McLure 
& Dunlop| ( |2004] l. 



2 Our Fell template from Vestergaard & Wilkes 1 2001} includes a 
small fraction of the Balmer continuum from I Zw 1. This is not the 
case of the template constructed by Tsuzuki et al. ( 2006 1 used by Wang 
eraL| ( |2009) to derive Eq.[T] 
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Fig. 2. Example of fit of the Mgll line for SDSS J0246-0825 
(black solid line) with a model (red) using the procedure de- 
scribed in Sect. |3.1| The model is the sum of a power-law con- 
tinuum (blue), a Fe II template (green), a broad Mg II component 
(dotted gray) and a narrow Mg II (dashed gray). The gray crosses 
indicate the region of the spectrum affected by atmospheric ab- 
sorption and not used for the fit. The dashed black line shows the 
residual of the fit (data-model). 



3.2. Lens modeling and intrinsic luminosities 

In order to derive the macro-magnification of the lensed quasar 
images of our main sample, we have modeled the lensing galaxy 
with an isothermal mass distribution using lensmodel vl.99o 
( Keeton|2001 1. The doubly-imaged systems have been modeled 
with a Singular Isothermal Sphere (SIS) and the quads with a 
Singular Isothermal Ellipsoid (SIE). For both doubles and quads, 
we accounted for the lens environment with an external shear 



term y. We used the same modeling technique as in Chantry 
et al. ( 20 1 0) > with slightly different observational constraints for 
the quads. In particular, we constrained the model with the rel- 
ative astrometry of the lensed images and lensing galaxy, and 
with the flux ratio of the lensed images. We used the flux ratios 
M derived from our spectra (see Sect.|4]and Table[5]l, accounting 
for 10% error bars on the flux ratio. Contrary to what we did in 
Chantry et al.| ( |20T0| , the flux ratios from our spectra have also 
been used to model the quads. We provide in Table [2] the values 
of the convergence k, shear y, shear position angle 6, and macro- 
magnification M (with the sign reflecting the image parity) pre- 
dicted by the models at the location of the lensed images. Note 
that for two objects, published HST-based astrometry was lack- 
ing and we used instead CASTLES^] measurements. We should 
mention that the flux ratios derived for the doubly imaged quasar 
SDSS 0246-0825 cannot be reproduced with a SI S+y model. 
This result was also found in the discovery paper (|Inada et~aL 



2005 ) and suggests that a companion object may strongly mod- 



ify the lens potential. The visible companion galaxy Gl located 
~ 1.8" from the lens cannot produce the anomaly and another in- 
visible companion has to be included in the model. For this sys- 
tem, we searched for models with a second SIS located +1" from 



the position of the putative companion derived by |Inada et al. 
fl2003] l. 



|http : //www . cf a . harvard . edu/castles/ 
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References: (1) Chantry et al. (2010), (2) Lehar et al. 1 2000 1, (3) Ina da et al. 
j2005) , error bars increased to 0.003"due to saturation in the PSF, (4) |Courbrn 
eTaLlpOTl), (5)|Sluse et al.H2008b) , (6) |Sluse et a^pm) , (7) CASTLES, (8) 
Vuissoz et aLH2008t 

Table 2. Lens models for each system of our main sample: 
convergence k, shear y, shear position angle 9 and macro- 
magnification M at the position of the lensed images. The ref- 
erences correspond to the lens astrometry used to perform the 
model. 



For the extended sample (Sect. |2.3[ ), we have not calcu 
lated new lens models. Instead we used the following macro 



magnifications factors M J A WM - 26.2 (Keeton et al. 2006), 



et al. 



2009jl, Mf 1 ^ 



3.9 (Kochanek 
We converted intrinsic 
nosities using Lboi 



2004) 



M j im = 11.6 ( Sluse et al.||2012|, = 10.3 ([MacLeod 



luminosities into bolometric lumi- 
BC x Lref, where BC is the bolomet- 
ric correction and L re f is the reference wavelength. We used 
the bolometric correction BC=(3.81, 5.15, 9.6) corresponding 
to Lref =(Li35o, £3000,^5100) ^ rom Sh en et al] ( [200 8 ). We used 
L re f = L3000 when we detect Mgll and Z^f = L5100 when we 
observe HB. We used L le f = L1350 in the other cases. 

3.3. Size of the BLR and Einstein radius 
We calculate the size /?blr of the BLR using the virial relation: 
GM BH 



R 



fFWHM 2 



(2) 



where G is the universal constant of gravitation and / is a fac- 
tor of order of one which encodes assumptions regarding the 



BLR geometry ( Vesterg aard & Peterson|2006||Collin et al.|2006 
|Decarli et al.|2008||La Mura et al.|2009|]Graham et al.|201 1) . As 
we perform only relative comparison between objects, the exact 
value of / does not matter. We therefore assume the same value 
/ = 1 for the whole sample. Because occurrence of ML does not 



only depend on the source size but also on the Einstein radius of 
the lens, we calculate Rblr/Re, where Re is the angular Einstein 
radius of a microlens projected onto the source plane as 



Re = 



AG (M) D 0S D h 



D„ 



(3) 



where the D's are angular diameter distances, and the indices o, 
/, s refer to observer, lens and source. (M) is the average mass of 
microlenses that we assume (M) = 1 Mq. 

3.4. Redshifts 

For nine out of thirteen objects re-analyzed here, we were able 
to measure the redshift of the source based on the Mg II emis- 
sion line, which is thought to give a good proxy towards the 
systemic redshift of the object. We found measurement com- 
patible with literature data (i.e. within 5z = 0.003) for six ob- 
jects (SDSS J0806+2006 and SDSS J1335+0118, SDSS J0246- 
0825, SDSS J1226-006, Q1355-2257, WFI2033-4723). The dif- 
ferences for the other objects probably arise from the use of dif- 
ferent emission lines used for the redshift calibration. This is 
clearly the case of HE 0435-1223 whose redshift was measured 
based on the C IV line which is known to be prone to systematic 
blueshifts in many quasars. The redshifts we measured based on 
Mg II are reported in Table [T] 

3.5. Physical properties of the sample: final remarks and 
uncertainties 

We present in Table [3] the physical properties of all the lensed 
sources discussed in this paper. Because we only measure differ- 
ential ML between two macro-images, we report values of Lboi.t 
and Lboi,2 derived for each of the lensed image. We use the aver- 
age between L\ and L2 for the calculation of Mbh and L/LEdd- 

The errors on the quantities reported in Table [3] were cal- 
culated in the following way. The error on the luminosity L is 
caused by a) microlensing, b) uncertainty on the absolute image 
magnification associated to the lens model, c) flux calibration 
and, d) intrinsic variability. We assumed that the microlensing 
budget error can be estimated from the spread between L\ and 
L2 and that the other sources of errors (b+c+d) correspond to 
20% of the flux of L. We added these errors quadratically. The 
calculation of the error on the FWHM is explained in Sect. |3.1| 
In case of measurements from literature, when no error bars on 
the FWHM were provided, we conservatively assumed cr ewhm 
= 0.25 FWHM. The errors on the other quantities were derived 
using error propagation formulae. When we report M BH from 
literature, we used a typical 0.4 dex error when no error bar was 
reported. 



4. Microlensing analysis of the spectra 

Several techniques can be used to unveil microlensing (ML) in 
lensed quasars. The most common method is the measurement of 
the emission lines equivalent width but this mostly allows one to 
assess microlensing of the continuum. Alternative methods ex- 



ists such as "multi-variable" scaling of the spectra (e.g. Angonin 
|et al.|1990| |Burud et al.||2002| |Wucknitz et al.|2003) , the tech- 
nique introduced by Popovic & Chartas ( 2005 ) based on multi- 
epoch spectra, flux ratio measurements after a multi-component 
decomposition of the spectra ( |Wucknitz et al.|2003[ |Sluse et aL 
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Table 3. Physical quantities associated to the source quasar. Values in parentheses are from Peng et al. 
( 201 l| l. The last four systems under the horizontal line are part of the extended sample which gather data from literature (Sect. 2.3 i. 



2007 1 Eigenb roTetlT1[2008]l or the so-called Macro -micro de- 
composition ( |Sluse et al. 2007| l. The latter technique, that we 



recall in Sect. 4.1 is the best suited for our goals because it en- 
ables an easy visualisation of the differential ML affecting the 
QSO spectrum, including partial ML of an emission line. It also 
provides a reliable estimate of the amount of ML affecting the 
continuum and of the intrinsic flux ratio M between the lensed 



images (see |Sluse et al.|2007| |2008a| |201 1| |Hutsemekers et al. 



2010 for comparison with other methods). 



4.1. Macro-micro decomposition (MmD) 



Following |Sluse et al] ( |2007[ > and [Hutsemeker s et al.| ( |2010| l, 
we assume that the spectrum of an observed lensed image F, 
is the superposition of 2 components, one component F M which 
is only macro-lensed and another one, Fm^, both macro- and 
micro-lensed. According to this procedure, it is possible to ex- 
tract F M and F Mfl by using pairs of observed spectra. Defining 
M — Mi IM2 (> 0) as the macro-magnification ratio between 
image 1 and image 2 and ju as the relative micro-lensing factor 
between image 1 and 2, we have: 



F\ = M x F M + M x // x F M fi 
F2 = Fm + Fmh ■ 



(4) 



To extract Fm and Fmh when M is not known a priori, these 
equations can be conveniently rewritten 



Fm - jzm{^-F 2 ) 

Fmh - J^m ( M ~ Fl ) ' 



(5) 



where A — M x fi. As explained in more detail in Hutsemekers 
et al. (2010, Sect. 4.1 and Appendix A), the factor A can be ac- 
curately determined as the value for which Fm(A) = in the 
continuum adjacent to the emission line. The factor M is cho- 
sen so that there is no visible emission above the continuum in 
Fmm(M) at the wavelength of the narrow emission lines which 



originate from regions too large to be microlensed. Often, nar- 
row emission lines are absent of the observed spectrum and then 
M is chosen to minimise the emission above the continuum in 
Fm m at the location of the broad line (see Sect. 4.2 for an il- 
lustration). If only a portion of the broad emission line is not 
microlensed, then M can be determined using this portion of the 
line (see Hutsemekers et al.|2 010). The micro-(de)magnification 
factor /j is simply derived using fi = A/M. 

The factor M is generally wavelength dependent because dif- 
ferential extinction between the lensed images produced by the 
lensing galaxy may take place. In addition, chromatic ML some- 
times occurs due to wavelength dependence of the source size, 
the blue continuum emitting region being smaller and more mi- 
crolensed than the red one, we may also have a wavelength de- 
pendence of //. For this reason it is necessary to estimate A and 
M in the vicinity of each emission line. Because differential ex- 
tinction affects both the emission line and the continuum, it does 
not modify the decomposition outlined in Eq.|5]as far as A and 
M are derived at the location of each line. 

Our method should ideally be applied to pairs of spectra sep- 
arated by the time delay. In Appendix [A] we show that intrinsic 
variability generally leads to small errors on M and fi and could 
rarely mimic microlensing deformation of the emission lines, 
only when the time delay between the lensed images is large 
(typically >50 days). Consequently, we apply the MmD to spec- 
tra obtained at one epoch, like those described in Sect. |2j with 
the caveat that in rare occasions intrinsic variability could be an 
issue. 



4.2. Example 

We illustrate in this section how the MmD is used to derive A, M 
and fi and to unveil microlensing of the emission lines. In Fig. [3] 
we show the MmD applied, at the location of the Mg II emis- 
sion, to the spectra of images B and D of HE 0435-1223. The 
MmD is performed in the following way. Up to a scale factor, 
Fm only depends on A, and Fm/i on M. Therefore we first mea- 
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Fig. 3. Macro-micro decomposition (MmD) applied to the 
Mg II line in HE 0435-1223 using the spectra of images B (Fi) 
and D (F2). The bottom pannel shows Fm for three different val- 
ues of A and an arbitrary value of M. The best value of A is 
A ~1.2 because it leads to Fm = in the continuum regions 
blueward and redward of the emission line. The upper pannel 
shows the decomposition for 3 different values of M. The best 
value is M =1.47 because it minimizes the emission in the line, 
i.e. keeping the flux at the location of the line just above the ap- 
parent local continuum depicted as a dotted black line. 



sure the factor A by making F M — in regions of the continuum 
blueward and redward of important emission lines, in this case 
Mg II (cf. Eq.[5]l. Which spectrum corresponds to F\ or F2 is ar- 
bitrary. We choose F\ = F B and F 2 = F D . In the bottom pannel 
of Fig. [3] we see that A ~ 1.2 leads to Fm ~ in the continuum 
regions. Alternatively, it is also possible to measure A in the con- 
tinuum regions from F1/F2. If the values blueward and redward 
of the line are different, we average them. Second, the factor M 
is chosen so as to minimize the emission above the continuum 
in Fmu at the location of the broad lines. The upper pannel of 
Fig.|3]shows the decomposition for three values of M. The value 
M — 1.37 is not adequate because it leads to a dip in Fj^ in the 
range 2750-2800A. Only values M > 1.47 lead to a valid de- 
composition. Finally Fm and Fj^ M are scaled according to Eq. [5] 
so as to have F2 = Fm + Fm^- The amount of microlensing // is 
derived with fi = A/M. For the example of Fig. [3] we find (A, 
M, yu)=(1.2, 1.47, 0.82). Larger amplitudes of ML are conceiv- 
able provided M is modified such that A = M Xfi still holds. We 
also emphasize that, with only a pair of spectra, it is not possi- 
ble to know if image #1 is magnified (resp. de-magnified) by 
or if image #2 is de-magnified (resp. magnified) by 1/fi. The er- 
rors on A and M are determined from the range of values which 
provide acceptable solutions to Fm and Fu^ The quality of the 
spectral decomposition, and errors strongly depend on the S/N 
of the spectra and on the strength of the microlensing effect (i.e. 
if — > 1, then M — > A in Eq. 15). Another illustration of the 
method is proposed in Appendix [B|based on simulated spectra. 



4.3. Results 

The MmD at the location of the main emission lines is shown in 
Fig. [4] The spectra of the systems Q0 142- 100 and HE 2149-2745 
from our main sample are not shown because of the absence of 



microlensing. A detailed description of the microlensing signal 
is provided in Appendix|C (resp. Appendixp) for the main (resp. 
extended) sample. In Table |4j we summarize the results of the 
MmD. For each object we provide A, M, fi, we list the possible 
sources of chromatic variations of A and the part of the broad 
emission line affected by microlensing, i.e. the line core (LC) 
which corresponds to the velocity range [-500: +500] km/s, the 
blue wing (BW), or the red wing (RW). When the signal is uncer- 
tain (because it is weak or depends on the exact value of M), we 
put the value in parentheses. When we do not detect microlens- 
ing of the continuum, we do not fill these three columns. The 
definition of these three regions does not rely on a strict defini- 
tion based on the line width and velocity ranges because of the 
variety of lines studied (sometimes blended with other lines) and 
of the uncertainty of the systemic redshift. We found, however, 
that this qualitative rating gives a fairly good synthetic view of 
the observed signal, with the drawback that it does not reflect 
the relative amplitude of each microlensed component. A more 
quantitative description, based on multi-component line fit and 
simulations, will be presented in a forthcoming paper. 



4.4. Accuracy of the intrinsic flux ratios 

We check that our measurement of the flux ratio M is a good 
estimate of the macro-model flux ratio. For this purpose, we 
compare it with flux ratios collected from literature, measured in 
the near-infrared or at larger wavelengths. Near-infrared flux ra- 
tios are believed to be less affected by micro-lensing because of 
the larger size of the emitting region (Wambsganss & Paczynski 
|1991[ ). They should therefore be a good proxy of the macro- 
model flux ratio. However, we may not exclude that in some 
cases, the size of the emission region is small enough to be af- 
fected by microlensing ( |Agol et al.|2009||Kishimoto et al.|2009| 
|201 p . Table[5jsummarizes measurements of the flux ratios from 
our spectra (col.#3) and at larger wavelengths (col.#4, from lit- 
erature). In that table, the value of M in col. #3 is an average 
between the values derived in the blue and the red part of the 
spectrum when we cannot disentangle between chromatic ML 
and differential extinction (the spread is then quoted in paren- 
theses and becomes our formal uncertainty). Otherwise, the val- 
ues are measured at our reddest wavelength. Figure [5] shows the 
pseudo-color ARm = -2.5 log(M w /M), where M is the flux ra- 
tio derived from the MmD and M K a broad band measurement 
in various spectral window (Table [5]). Three spectral windows 
are considered. First the region around 7500 A (which corre- 
sponds to the reddest part of our spectra). Second, the //-band 
(~ 1.6 jUm). Third, the reddest available wavelength (ZT-band 
= 2.2 yum, L-band = 3.8 /mi, 11.5pm or radio 8.4 Ghz; see 
Table [5J. There are five objects which are lacking flux mea- 
surements at a wavelength redder than //-band. In that figure, 
we have quadratically added the expected variability over the 



time delay period (see Table C. 1 in Appendix) to the error bar of 



the continuum flux ratio. This error in principle affects M, and 
should therefore be propagated to all the other measurements. 
However, to ease legibility and identify more clearly discrepan- 
cies, we have included the variability error budget only to the 
spectroscopic flux ratios. 

Figure [5] shows that there is a general trend towards AR M = 
with increasing wavelength (black circle symbols), confirm- 
ing that our estimate of M is a good proxy to intrinsic flux ra- 
tios (see Appendix [C] & |D]for a discussion of each object in- 
dividually). Small discrepancies of up to 0.1 mag are however 
observed because of intrinsic variability, and possibly of dif- 
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Fig. 4. Macro-micro decomposition (MmD) around the main emission lines. In each pannel, we show the part of the line profile 
which is only macrolensed, Fm, and the part of the line profile which is both macro- and micro-lensed, Fm m . The spectrum of Image 
#2 (F2 = Fm + FMfi) is superimposed (dotted line). The vertical dotted line denotes the rest wavelength of the emission line (using 

Fm/i contains the micro-magnified continuum 
le emission lines unaffected by microlensing. 



z s from Table[T|. The values of M and fi are given for each decomposition (Sect. 4.1 
and, quite often, a part of the emission line profiles, while Fm contains the bulk of t 
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CML, (DE, IV) 
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(r) Q2237+0305+ 


A/D 
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CIV 


0.425+0.002 
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i.ou_ l3 


0.425+0.085 


(e) 
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0.438+0.003 


l 00° 2,i 
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1 



Notes: (a) Uncertain due to the proximity to the edge, (b) Likely cross contamination between the lens and the QSO images, (c) The noisy decomposition around 
CIII] might hide a faint broad component similar to the one observed in Mgll. There is a narrow component in the red wing superimposed to the broad one. (d) 
Spectra already corrected for chromatic effects induced by the host galaxy (see Sluse et al. 20071. (e) Spectra already corrected for differential extinction between A 
& D (seeiSluse et al. 201 111. t Objects from the extended sample. 

Table 4. Results of the MmD technique (see Sect. 4.3 and Appendix |c| & [o]for details). Cols. #1 and #2 give the object name 



and the image pair used for the decomposition, cols. #5, #6, and #7 give the continuum flux ratio A, the macro-magnification ratio 
M, and the micro-magnification /j. cont of the continuum. These quantities are given at the approximate wavelength (observed frame) 
given in col.#3, which corresponds to the center of the emission line given in col. #4 (see Sect. |4.1| >. Column #8 lists the possible 
origins of the chromatic changes of A, namely differential extinction (DE), chromatic microlensing (CML), contamination by the 
lens (LC), or intrinsic variability (IV). The last 3 columns indicate if microlensing of the broad line is seen in the blue wing (BW), 
line core (LC; [-500: +500] km/s), or red wing (RW). The value is used when there is no microlensing and 1 when it is present. 
The value is in parentheses when the signal is weak or depends strongly on the exact value of M. 



ferential extinction. The effect of the intrinsic variability is in 
general smaller than the prediction from the structure function 
(Table 



C.l 



Among the 12 systems for which K, L, MIR or ra- 
dio flux ratios were available, five objects (depicted with open 
symbols) have flux ratios in the red which disagree significantly 
with our estimate of M. In one case (J 1 131-1231), it is likely that 
the flux ratio in K-band is significantly affected by microlens- 
ing. Two other systems (HE 0435-1223 and SDSS J0806-2006) 
are already flagged as anomalous becaus e of the discrepancy be 



tween their K- and L-band flux ratio^](Fadely & Keeton 201 1 



|2012 1. This anomaly might be produced by a massive substruc- 
ture in the lens ( Fadely & Keeton||2012 1, and deserves further 
investigation. The origin of the disagreement between our esti- 
mate of M and the K-b&nd flux ratios for the last two objects 
(SDSS Jl 138+0314, SDSS J1335+0118) is less clear. The most 
simple explanation is intrinsic variability larger than our predic- 
tions and/or significant differential extinctiorr] 

Figure [5] also allows us to test whether ML and differential 
extinction are negligible in //-band, as often stated in literature. 
We clearly observe in the figure that the //-band value of ARm 



4 The case of HE 0435-1223 is a bit puzzling since two discrepant 
measurements of the A"-band flux ratios have been published by Fadely 



& Keeton|j2(5TT| and |Blackburne et al.|j20T7) . 



For SDSS Jl 138+0314, we cannot exclude a systematic error due 
to the slit centering because we could not cross check the spectroscopic 
flux ratio with broad band simultaneous data. 
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Fig. 5. Pseudo color ARm - -2.5 log(M„,/M), where M is the intrinsic flux ratio derived from our spectra using the MmD and M w 
is the flux ratio measured in 3 different wavelength ranges: (1) the region around 7500 A (blue diamonds), (2) the //-band (red 
triangles), and (3) wavelengths redder than //-band (black circle). Open symbols are used for the systems where our spectroscopic 
estimate of M disagrees with the value at wavelengths redder than //-band (see Sect. 4.4 1. The objects are ordered from left to right 
as a function of increasing microlensing of the continuum (note that the steps are not linear). The gray labels provide the name of 
the systems (shortened to the first 4 digits of the RA) and the upper axis indicates the redshift of the source. 



is systematically closer to (no ML) than the value at 7500 A. 
This trend supports the idea that intrinsic variability and differ- 
ential extinction do not add much noise in the estimation of M 
(otherwise the effect would not be systematic). In addition, ARm 
is in general less close to in //-band than at redder wavelength. 
This confirms that ML is often detectable in //-band. We would 
expect the source at larger redshift (i.e. for which //-band cor- 
responds to the bluest and therefore smallest emitting regions) to 
be more prone to ML in the //-band than the low redshift ones, 
but //-band ML is also observed in the low redshift sources. 
There is neither a clear correlation with the black-hole mass, 
which would be naively expected because the source size scales 
with the black hole mass. The absence of such correlations is 
likely related to ML variability, because the //-band data and 
optical spectra are not simultaneous, and to the variety of mi- 
crolensing events which may not all lead to significant ML in 
the //-band. 

The above discussion confirms the efficiency of the MmD to 
derive intrinsic flux ratios, with a typical inaccuracy of 0. 1 mag 
associated to intrinsic variability and, in a few cases, systematic 
errors associated to differential extinction and microlensing. 

5. Discussion 



the occurrence of microlensing of the broad emission lines and 
try to identify some trends as a function of the object physical 
and spectral properties (Sect. 5.2 1. We finally discuss the conse- 



The objects of our main sample (Sect. 2.1 1 were targeted with 
the goal of measuring the redshift of the lensing galaxy and in- 
dependently of any known detection of microlensing. Therefore, 
we can make some basic statistics regarding the chance of de- 
tection of microlensing in a lensed quasar based on these sys- 
tems. Because observations were not performed exactly at the 
same time but over a two-years period, and because of possi- 
ble secondary selection biases (e.g. selection of the brighter tar- 
gets, large source redshift range, bright lensing galaxy, mix of 
doubles and quads), this estimate should not be very accurate. 
On the other hand, we can use the complete sample to discuss 
the variety of microlensing deformations of the broad lines. We 
first address the occurrence of microlensing of the continuum 
and its chromatic variations in Sect. |5. 1 1 Afterwards, we discuss 



quences of our observations for our understanding of the BLR in 
Sect. 1331 



5. 1 . Microlensing of the continuum 

Our observations confirm that microlensing of the continuum 
is observed at a level > 0.05 mag in 11 out of 13 systems, i.e. 
~ 85% of the systems. Because Q0142-100 and HE2149-2745 
are the 2 systems with the largest black hole mass, we may be 
tempted to associate the absence of microlensing in these sys- 
tems to a large size of the continuum emitting region. However, 
simple estimates of the latter based on the accretion disc theory 
lead to continuum sizes at least 10 times smaller than the mi- 
crolens Einstein radius (Mosquer a & Kochanek|[20TT] l. On the 
other hand, the past detection of microlensing in HE21 49-2745 
( |Burud et al.| 2002) teaches us that microlensing of the contin- 
uum can be observed in these massive quasars. 

Despite the large occurrence of microlensing of the con- 
tinuum, significant chromatic variations over the optical wave- 
length range, and more specifically those possibly caused by mi- 
crolensing, are less common. Only four out of twelve systems 
show changes possibly associated to microlensing (HE 0047- 
1756, SDSS J0246-0825, HE 0435-1223, Q1355-2257). This 
probably reflects that the source often lies far enough from 
a microlensing-caustic such that chromatic ML in the op- 
tical range is only weak. We should however notice that 
microlensing-induced chromatic changes are expected for all the 
microlensed sources if the wavelength coverage is large enough 
(i.e. extending up to 2.5 /zm, cf. Bate et al. 2008 ; Yonehara et al. 
|2008l|FToyd~et al.|2009[|Blackburne et al.|201 1) . 



5.2. Microlensing of the BLR 

The most common broad emission lines (BEL) detected in our 
sample are C III] and Mg II. Therefore we focus our discussion 
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- The CIII] and Mgll emission lines are simultaneously 
observed in nine systems. For two of them, microlens- 
ing affects none of the lines and for six others, both lines 
are affected (but the signal in CIII] is uncertain in two 
cases). Microlensing of CIII] only, is seen in one system 
(WFI 2033-4723), but in this particular case the signal might 
also be due to microlensing of the Al III or Si III] lines 
blended with the blue wing of C III] . 

- In most cases either the red wing or the blue wing is mi- 
crolensed, with roughly equal occurrence. In addition, the 
microlensing signal always starts relatively close to the sys- 
temic redshif|^](typically within 1000 km/s). 

- Two of the six systems where microlensing of Mgll and 
CIII] is simultaneously observed display Fm/i with a 
different shape in Mgll and CIII]. These systems 
(SDSS J1335+0118 and Q1355-2257) are also those ones 
where a significant blueshift (~ 1000 km/s) of C III] is visi- 
ble. 

- For one system (SDSS 11335+0118), the microlensing sig- 
nal suggests two velocity components in the BLR, one which 
gives rise to a broad symmetric profile and the other one to a 
narrow component centered at v~ 1500 km/s. These compo- 
nents are detected in both C III] and Mg II profiles despite 
their different shapes. Note also that the Mg II line in this 
object is more asymmetric than in any other system of our 
sample. 

- The strongest microlensing effects in the BEL are in gen- 
eral associated to a symmetric signal affecting both positive 
(red) and negative (blue) velocities (Q2237+0305, 11131- 
1231, BRI0952-115 and SDSS 11138+0314). This effect is 
expected for most of the BLR geometries studied in the liter- 
ature when the source lies in a large demagnification valley 
of the micro-caustic network or when a large fraction of the 
source crosses a micro-caustic (e.g. Abajas et al.||2007} . In 
the first case the symmetry should last several crossing time- 
scale^] while in the second case, the effect should be more 
transient. 

- Figure [6] shows the properties of the microlensing signal in 
Mgll or CIII] as a function of the bolometric luminosity 
Lboi an d °f Rblr/Re- Almost any kind of microlensing sig- 
nal is observed at a given Lboi or Rblr/Re- The objects with 
smaller Rblr/Re seem more prone to stronger deformation 
of the lines, as expected. Possibly because of the large error 
bars, we did not detect any correlation of the ML signal with 
the other physical properties calculated in Table [3] 



(17)|Eigenbrod et al.|(|2006aj>r(18)|Sluse et al.|(|2006^; (19)|Falco et al.| 5.3. Consequences for the BLR 



( 1996 ); (20) |Chantry & Magaml ( |2007| l; (21) |MacLeod et al.H2009> 
Table 5. Comparison of M from this paper (col. #3) with values 
derived in the literature (col. #4). The object name (col. #1) has 
been shortened to the first 4 digits of the RA. The last column 
indicates the filter/line for which the literature measurement has 



been performed. See Sect. 4.4 for discussion. 



on these two lines. Considering only the objects showing mi- 
crolensing of the continuum in our main sample, we found mi- 
crolensing of one of these two lines in ~ 80% of the systems 
(8 out of 10). We now extend the discussion to the whole sam- 
ple and address the question of the variety of line deformations. 
Looking to the results of the MmD we found: 



The current analysis has identified various signatures of mi- 
crolensing of the BLR of quasars. A confrontation of the ob- 
served signal to different geometries of the BLR is beyond the 
scope of the present paper. However, the data presented here al- 
low us to derive several qualitative results. 

(a) The typical amount of microlensing we observe suggests 
that the BLR cannot be much more extended than the mi- 
crolensing Einstein radius Re- Assuming that 10% of the 
emission line flux is microlensed (see Fig. [4]), the volume of 
the region magnified by a microlens of R E should roughly 
constitute 10% of the total emitting volume. For a spherical 



This is not the case for J1226-0006 but the signal is flagged as pos- 
sibly contaminated by the flux from the lensing galaxy. 

7 The crossing time scale is the time it takes the source to cross a 
distance equivalent to its own radius 
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BLR, this means nR\ x 2R BhR ~ 0.1 x 4/3 n x R^ lr , i.e. 
Rblr ~ 4/?£. A similar value is obtained for a face-on disk. 
Although this is a very simplified treatment of the problem, 
this is in agreement with the independent estimates of Rblr 
given in Table [3] 
(b) The frequent observation of asymmetric blue/red deforma- 
tions of the emission lines demonstrates that the BLR does 
not have, in general, a spherically symmetric geometry and 



velocity field. As shown by the simulation of Schneider & 
|Wambsgans"sl ( |1990| l; |Abajas et ah] ( |20Q2) >; |Lewis & Ibata 



(2004), asymmetric deformation of the broad fines can only 



be produced when the BLR has an anisotropic geome- 
try or velocity field (i.e. axially symmetric, biconic, ...). 
Note that the symmetric microlensing of the broad line ob- 
served in some systems (e.g. Q2237+0305, Jl 131-1231, 
SDSS 1 138+0314) may still occur with an anisotropic BLR. 
In this case the signal could become asymmetric and re- 
vealed by a dedicated monitoring campaign. 

(c) The comparison of the microlensing signal in Mgll and 
C III] suggests that differences observed in emission lines 
of different ionisation degrees are more pronounced when 
one of the two lines is not centered at the systemic red- 
shift. This may be of particular importance to disentangle 
blueshift caused by obscuration of the redward part of the 
line and blueshift associated to the dynamics of the BLR 
(e.g. |Richards et al.|2002|[20TT) . 

(d) For one system (SDSS J0246-0825), a symmetric microlens- 
ing of the emission line is observed with a dip in the line 
core. This is strikingly similar to the signal expected for a bi- 



et al. 



Micro" 



conical disk ( Abajas et al. 2002[ |Lewis & Iba ta 2004 Abajas 



2007l. However, the solution may not be unique. 



ensing of keplerian disks (modified disks like those 



vational bias (microlensing of the other wing is too small 
to be detected), or not. In this case, it seems, according to 
published simulations, that microlensing of a keplerian disc 
qualitatively better reproduces the observed signal. 



6. Conclusions 

In this paper we have searched for the presence of microlens- 
ing, in the continuum and in the broad emission lines, among 
strongly lensed quasars based on optical spectroscopy of a sam- 
ple of 13 systems (3 quadruply and 10 doubly-imaged systems). 
The spectra of these systems, originally targeted for detecting 
and measuring the redshift of their lensing galaxy (Eigenbrod 
|et al.|[2006b| [2007), have been re-reduced and deconvolved in 
order to accurately deblend the flux of the lensed images and 
of the lensing galaxy. For the three systems BRI 0952-0115, 
SDSS J1138+0314, SDSS J1226-0006, these are the first pub- 
lished spectra of the individual lensed images. In order to get a 
more complete overview of the variety of microlensing signals 
which can be detected, we have complemented our main sam- 
ple with previously analysed objects, i.e. SDSS J0924+0219, 
HE 2149-2745, Jl 131-1231, Q2237+0305 and H1413+117, ex- 
tending our sample to 17 objects. 

We derive robust estimates of the intrinsic flux ratios M be- 
tween the lensed images and of the amplitude fi of microlensing 
of the continuum. Based on optical spectra only, we have shown 
that we can retrieve M to typically 0.1 mag accuracy. Higher 
accuracy may be expected by using spectra separated by the 
time delay, over a wavelength range extending from UV to NIR. 
Our ability to derive accurate M and fj. is particularly relevant 
for the study of the so called "flux ratio anomalies" produced 



introduced in Abaj as et aL|2002| or relativistic discs as dis- by dark matter substructures or microlensing ( |Mao & Schneider 



cussed in Popovic et al. 2001 1 sometimes lead to a dip in 



the microlensed signal at zero velocity. Only one other sys- 
tem (HE 0047-1756) shows ML similar to the one observed 
in SDSS J0246-0825 but with a larger signal in one wing 
than in the other, a case which should in principle be more 
likely (Abajas et al. 2002 2007 1. This might be an obser- 



1998; Dalai & Kochanek 2002). Indeed, current studies of this 
effect are challenged by the limited number of reliable flux ra- 
tios HKeeton et al.|2005l|Metcalf & Amara|20T2"l|Xu et al.|2012| i. 
Therefore, the use of our technique should help in increasing the 
number of interesting systems. Our results may also be used to 
derive the fraction of smooth matter in lensing galaxies. A study 
of the probability of deriving the observed values of yu in sample 
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of lensed quasars, varying the fraction of compact/smooth dark 
matter along the line-of-sight of the lensed images, makes pos- 
sible to constrain the fraction of matter in compact form in these 
galaxies (|Schechter & W ambsganss 2 004| |Pooley et aL [2009 ; 
Mediavilla et al.|2009||Bate et al.|201 l||Pooley et al.|2012| i. 



We find in our sample, which contains sources with bolomet- 
ric luminosities in the range [10 44 - 7 , 10 47 - 4 ] erg/s and black hole 
masses in the range [10 7 6 , 10 9 8 ] Mq, that 85% of the sources 
show microlensing of the continuum. The microlensing of the 
continuum is not systematically associated to significant chro- 
matic changes in the optical range. These chromatic changes are 
not necessarily absent but relatively weak. Because we observe 
microlensing of the continuum in our sample in the range 0.2- 
0.8 mag in /?-band (observed frame), microlensing may in gen- 
eral not be neglected in the //-band (1.6 fim). This implies that 
studies of the accretion disc based on microlensing will strongly 
benefit of a large wavelength coverage, extending at least up to 
the Zf-band. Another consequence is that microlensing may of- 
ten exceed differential extinction in the NIR, leading to possi- 
ble biases in studies of extragalactic extinction curves in lensed 
quasars. 

Using our MmD decomposition technique, we have been 
able to unveil microlensing-induced deformation of the emis- 
sion lines independently of a detailed modeling of the quasar 
spectrum. Among the systems with a microlensed continuum, 
80% show deformations of at least one broad emission line. The 
two major characteristics of the signal are its relatively low am- 
plitude (typically 10 % of the line is affected) and its variety. We 
searched for correlation between the observed microlensing sig- 
nal and the luminosity, black hole mass and Eddington ratio of 
the microlensed quasars, but we did not find any clear correla- 
tion. Contrary to most previous observations, we frequently de- 
tect microlensing of either the blue or the red wing of the broad 
emission lines instead of a signal affecting roughly symmetri- 
cally both components. This simple observation implies that the 
BLR does not have a spherically symmetric geometry and ve- 
locity field, at least in those objects. Microlensing of only one 
wing of the line has been observed recurrently in the past for 
SPSS J1004+4112 (|Rlchards et al.|[2004| |Lamer et al.||2006| 
Gomez-Alvarez et al.|2006| . |Abajaset al.| ( |2007] > has shown that 
the signal observed in that system is compatible with microlens- 
ing of a biconical wind but that in general the two wings should 
be affected by microlensing, with different amplitudes. Such a 
signal is detected in only one system of our sample. This sug- 
gests that the simple biconical model is not adequate for the BLR 
of the lines we studied (CIII] and Mgll). A modified keplerian 
disc (Abaj as et al.|2002| seems to be a promising alternative as 
a "generic" BLR model. 
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Appendix A: Effect of intrinsic variability on the 
MmD 



The MmD technique described in Sect. 4. 1 should ideally be ap- 
plied to pairs of spectra separated by the time-delay in order cor- 
rect for any effect introduced by intrinsic variability. Hereafter, 
we investigate how intrinsic variability may affect the MmD 
technique applied to spectra obtained at a single epoch. For this 
purpose, we have decided to create mock lensed systems based 
on existing spectra of quasars observed at several epochs. The 
principle of our simulation is to pick pairs of spectra of an ob- 
ject at two different epochs t\, ti. A pair of spectra simulates the 
single-epoch spectra of two images of a macro-lensed quasar 
with a time delay At = t\ — ?2 and M =1. Then, we amplify 
the continuum of one image to simulate microlensing and apply 
the MmD. Since only the continuum is microlensed, we do not 
expect emission lines or part of them in Fm^, except possible 
contamination due to intrinsic variability. 

Specifically, we proceeded as follows. First, we used pub- 
licly available reverberation mapping data^ of Palomar-Green 
quasars. In this database, we choose pairs of spectra of the same 
object separated by a delay At in the ranges (a) 1-20 days, (2) 
20-40 days, (3) 40-60 days, (4) 60-100 days. Second, we artifi- 
cially microlensed the continuum of the first spectrum by a factor 
/j =1.5. Third, we applied the MmD, estimating automatically 
A and choosing M to minimize the flux in Fmh at the position 
of the H /3 lineJ^J Note that we had to restrict ourselves to pairs 
of spectra obtained with the same instrumental setup to avoid 
spurious line deformation introduced by variable spectral reso- 
lution. We show in Fig. A.l the result of this procedure at the 
position of the H/3 line for PG0052 (R BLR ~ 134 light days), 
PG0953 (7?blr ~ 151 light days), PG1613 (R BL r ~ 39 light 
days), PG0026 (7?blr ~ 1 13 light days). We also report the mea- 
sured fractional variation of the continuum (e c ) and of the line 
(e[) during At. This figure illustrates that in general the defor- 
mations of the emission lines caused by intrinsic variability are 
too weak to mimic microlensing and introduce a significant sig- 
nal above the continuum in Fm/j at the location of H/3. When 
the delay becomes large (typically > 40 days), it happens that a 
weak signal is detected in the emission lines (e.g. pannel (c) and 
(d) for PG0052, pannel (c) for PG0026). This happens when e, is 
large and when it differs significantly from e c . From this figure, it 
seems that differences of e; and e c by more than 10% are needed 
to introduce noticeable line deformations in the decomposition. 
In order to derive how frequent this situation appears, we have 
calculated e//e c as a function of At for the objects of the sample 
of |Kaspi et al. ( 2000 1. We did not find a clear change of e;/e e nor 
of the standard deviation cr e with the size of the BLR. Therefore, 
we report in Table |A. lj the average value of e;/e c together with the 
average value of cr e . We see in Table A.l that e//e c is in average 



equal to 1 with a scatter < 10% on periods corresponding to a 
time delay < 50 days. Since most of our targets have At < 50 
days (13 out of 17 targets), we may safely conclude that statis- 
tically, intrinsic variability is unlikely to mimic microlensing of 
the broad lines for such time-delays. The situation might be less 
favourable for the objects with At > 50 days, but only two out of 
four of these systems (Q1355-2257 and WFI 2033-4723) show 
possible microlensing of the emission lines. 



http : / /wi se - obs . tau . ac . i 1/ - shai/PG/ 
9 For those spectra, we are able to derive M directly from 
[O III] which does not vary, but we did not choose this option to fol- 
low the same procedure as in Sect. 4.1 
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Table A.l. Average value of fractional variation of the H/3 line 
and of the continuum (e//e c ) on periods Af (col. 1), for the sam- 
ple of reverberation mapped quasars of published in Kaspi et al. 
(2000). The average scatter cr € is given in col. 3. 



Appendix B: MmD applied to a simulated spectrum 

Similarly to the example of MmD applied to HE 0435-1223 in 
Sect. 4.2 we show in Fig. B.l the MmD applied to mock spec- 



tra roughly mimicking our spectra of HE 0435-1223. The mock 
spectra of HE 0435-1223 are defined in the following way: 

M D {hd F c + n' D E a + E h ) 
i (fi B F c +fi' B E a + E b ) 



^■mock 

F§ ock = Mi, « 



(B.l) 



where F c is the continuum emission, (Eb, E a ) are gaussian emis- 
sion profiles centered on (Ae,,, AeJ =(2798, 2803) A and with 
(FWHM Eb ,FWHM Ea ) = (5700, 2500) km/s. In this equation, 
the macro model magnification M and the micro-magnification 
H (for the continuum) and fj 1 (for the line) of individual images 
have been written explicitly. These quantities have been cho- 
sen arbitrarily such that M - M B /M D = 6.46/4.39 = 1.47 
(matching the macro-model), fi - HbIhd = 0.74/0.9 = 0.82, 
and // = yUg/yU^, = 0.85/1.0 = 0.85. We used a different micro- 
magnification factor in the line and in the continuum to account 
for the fact that lines are emitted in a region larger than the con- 
tinuum (i.e. hi closer to 1 than /i c ). We have also added a fake 
atmospheric absorption to the spectra of B and D in order to 
increase the similarity with the observed spectrum of H 0435- 
1223. Despite the similarity with the data, the model of the 
emission line of Eq. B.l should be considered only for illus- 
tration purpose. The MmD applied to these spectra is the same 
as the one discussed in Sect. |4.2| In absence of noise, the value 
M = 1 .47 and /U = 0.82 are retrieved. The component Fm^ shows 
that a red fraction of the emission line (corresponding to our 
input component E a ) is retrieved as microlensed, in agreement 
with our input model. 

Appendix C: Characteristics of the main sample 

We provide hereafter detailed notes about the characteristics of 
the spectra of each object of the main sample. For each target we 
discuss i) the quality of the deconvolution, which might intro- 
duce spurious signal, ii) the chromatic changes observed in the 
spectra, iii) the microlensing-induced deformation of the emis- 
sion lines, when a deformation is observed, and iv) important re- 
sults from literature which shed light on the source of chromatic- 
ity and confirm/infirm our flux ratio measurements. We provide 
in Fig. C.l the spectral flux ratios between the pairs of images 



and we discuss in the text four origins for the chromatic changes 
in the spectral ratios: differential extinction (DE), chromatic mi- 
crolensing (CML), contamination by the lens (LC), or intrinsic 
variability (IV). In order to quantify the possible systematic er- 
ror on our estimate of M caused by intrinsic variability, we also 
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Fig.A.l. Application of the MmD method for H{3 & [OIII] for mock lensed quasars simulated based on spectro-photometric 
monitoring data of Palomar-Green quasars. For each pair of spectra separated by a delay At, the continuum of the one of the spectra 
has been artificially microlensed by ju =1.5 and M =1 has been assumed. For legibility, only the fraction F Mfl (red solid line) and 
the reference spectrum (blue dotted line) are shown. Each row corresponds to a different object. The time-delay between the pairs 
of spectra is increasing from pannel (a) to (d). For each pannel, we provide M retrieved with the decomposition, the variability in 
the continuum e c and in the line e/. 



provide in Table C.l the amplitude of variation Am of an object 
over the time-spent of the time-delay. This quantity has been de- 
rived using the g-band structure functions (divided in 6 bins of 
M BH and L, cf. their Eq. 2 and Table 3) of |Wilhite et aT] ( [2008] l 
and the black hole masses and luminosities of TablefS] In the fol- 



lowing, the references to the NIR flux ratios from literature are 
not systematically given. They can be found in Table [5] 

(a) HE0047-1756: Deconvolution: The deconvolution is 
slightly less good than for other systems. The total flux left in 
the residual under the QSO images never exceeds 0.1% of the 
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C.l. Flux ratio for all the image pairs of our main sample. 
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Fig. C.l. continued 



flux of the QSO image. Although this is a very small amount of 
flux, this appears as a systematic feature suggesting that the PSF 
is less representative of the QSO images than in other systems. 
Chromaticity: The ratio B/A shows a monotonic decrease with 
increasing wavelength. This chromatic change is not due to con- 
tamination by the flux of the lensing galaxy. Indeed, in order 
to reconcile the shape of the spectrum in A & B, one has to in- 
voke a contamination of A (the brightest image) by ~ 4 times the 
measured flux of G and a contamination 10 times smaller for im- 
age B. Therefore, the two most likely explanations are DE with 
(M(blue), M(red), /j) = (0.231, 0.220, 1.173) and CML with (M, 
//(blue), Mred)) = (0.220, 1.260, 1.170). 

Broad Lines: Whatever the origin of the chromatic changes, the 



blue wing of Mgll is microlensed. Microlensing of the blue 
wing of C III] is also tentatively observed but this is more un- 
certain owing to the proximity of this line from the red edge of 
the spectrum. 

Notes: The observed monotonic decrease of B/A was also re- 
ported in the discovery spectra of Wisotzki et al. (2004). For 



these spectra, obtained in Dec. 2001 and Sept. 2002, we mea- 
sure A = 0.347 + 0.002 and M = 0.24 + 0.02 at the wavelength 
of C III]. The different value of A found in our spectra is caused 
by time-variable microlensing. On the other hand, M is compat- 
ible with our measurement, confirming that intrinsic variability 
may not be too large. The //-band flux ratio from |Chantry et alT] 
(2010), does not follow the chromatic variation of M observed 
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Fig. B.l. Macro-micro decomposition (MmD) applied to sim- 
ulated spectra of HE 0435-1223 mimicking the observed data. 
The decomposition is similar to the one showed in Fig. [3] The 
bottom pannel shows F M for three different values of A and an 
arbitrary value of M. The best value of A is A ~1.2 because it 
leads to Fm — in the continuum regions blueward and redward 
of the emission line. The upper pannel shows the decomposition 
for 3 different values of M. The best value is M =1.47 because 
it minimizes the emission in the line, keeping the flux above the 
apparent local continuum depicted as a dotted black line. 



in the spectra and is found 0. 1 mag smaller than the value de- 
rived at the level of C III] . This may be explained by intrinsic 
variability (Am ~ 0.07 mag) and/or microlensing of the H-band 
continuum. 

(b) Q0142-100: Deconvolution: Because of the inaccurate 
PSF, residual signal above the noise is visible in the deconvolved 
image. However, this flux amounts only to 0.02-0.2% of the total 
flux in the lensed images. The flux of the lens galaxy G, located 
only ~ 0.4" from B, reaches 20% of the flux in image B in the 
red part of the spectrum. This could lead to significant contami- 
nation of image B by the lensing galaxy G. 

Chromaticity: There is no differential ML between the contin- 
uum and the emission lines however the two spectra do not su- 
perimpose once scaled with the same magnification factor. This 
chromatic effect is probably caused by residual contamination 
of image B by flux from the lens since only ~ 13% of the ob- 
served flux of G is needed to explain the observed chromatic 
trend. Intrinsic variability could also play a role. We discard DE 
because it involves a larger reddening of image A which is far- 
ther away from the galaxy than image B. 

Notes: Our flux ratio is in good agreement with those obtained 
by Fadely & Keeton (201 1) in the K— and L'— bands. Koptelova 
|et al.| ( |2010) measured flux changes by 0. 1 mag over a period of 
100 days, which probably explain the flux ratio differences with 
//-band (Leh ar et al.|2000) . Color differences associated to IV 
might also play a role and have been reported for this system 
by |Koptelova et aL ( 2010 1. The stud y of the extinction in the 
lensed images by Ostman et al. ( 2008 ) disfavours significant ex- 
tinction in this system, in disagreement with |Falco et al.| ( fl999| l; 
lEliasdottir et al.| ( |2006) . 

(c) SDSS J0246-0825: Deconvolution: The deconvolution is 
very good. The small systematic residual detectable in the vicin- 
ity of the brightest lensed image A contributes to less than 0.05% 
of the flux of A. 



Chromaticity: The observed chromatic trend in B/A is compat- 
ible with DE (M(blue)=0.32, M(red)=0.34, fj, =0.76) or with 
CML (M=0.34, //(blue)=0.73, //(red)=0.76). Contamination of 
image A by the host galaxy, unveiled as a ring feature close to A 
in Inada et al. (2005) is plausible but is probably very low due to 



the slit clipping and to the relative faintness of this feature. 
Broad Lines: Whatever the origin of the chromatic changes, the 
broad component of CIII] is microlensed but ML of Mgllis 
more tentative. 

Notes: There is a good agreement between our spectral-based 
estimate of M and the H— and L— band flux ratios ( Inada ^eFaL] 
2005 ; Fadely & Keeton 201 1 ). The agreement in the /T-band is 
slightly less good but still marginally consistent with the other 
measurements (Fa dely & Keeton|201 l| l. 

(d) HE 0435-1223: Deconvolution: Very good results are 
obtained with the deconvolution. The lensing galaxy is relatively 
bright compared to the lensed images (only 4 times fainter than 
the flux of the lensed images above 6500A) and we may not 
exclude contamination of the latter by the lens. The symmetric 
location of the lensed images aside the lensing galaxy and their 
similar brightnesses however argue against large differential ef- 
fects (i.e. if contamination takes place, the spectra of A & D 
should be corrupted the same way by the lens). 
Chromaticity: We observe a chromatic increase of B/D with in- 
creasing wavelength, in agreement with the chromatic changes 
observed by [Fadely & Keeton ( 201 1 \ based on HST images. This 



trend may not be produced by contamination from the lensing 
galaxy because one needs a large contamination of image B (by 
at least 30% of the observed flux of G) and nearly no contami- 
nation of image D to mimic this effect. Because we find approx- 
imately the same microlensing factor in the blue and in the red, 
DE (M(blue)~1.34, M(red)~1.47, n -0.81) is the most natural 
explanation, B being more reddened than D by the dust in the 
lensing galaxy. Alternatively, CML may be at work (M ~1.47, 
//(blue) ~ 0.74, //(red) ~ 0.82). 

Broad Lines: Whatever the origin of the chromatic changes, ML 
of the red wing of the C III] and Mg II lines is observed. 
Notes: Two fairly different flux ratios in A"-band have 
been reported in literature. Fadely & Keeton (2011! reported 
B/D =1.49+0.12 while |Blackburne et al.| p0TT| reported 
B/D =1.27±0.04. The latter estimate agrees well with the 
//-band and L-band flux ratios while the former one agrees 
with our spectroscopic estimate M =1.47. Fadely & Keeton 



(201 1 2012| l interpreted the ZT-band flux ratio as caused by mi 
crolensing. However, the simplest microlensing scenario is hard 
to reconcile with the non monotonic variations of B/D, if real. 
Additional data are needed to solve this puzzle. 

(e) SDSS J0806-2006: Deconvolution: The deconvolution is 
good but the spectra have a significantly lower signal to noise 
than for the other systems. The residual flux under the point-like 
images reaches up to 0.2% of the flux of the QSO image. The 
lensing galaxy is well deblended from the QSO images. 
Chromaticity: The estimate of A at wavelengths shorter than 
CIII] is uncertain due to the proximity of the noisy edge of 
the spectrum. Despite of this, the measurements are consistent 
with a flat ratio B/A in the continuum from 4000 to 8000 A and 
ML affecting only the continuum emission. 
Notes: Our average macro-magnification ratio M B a ~ 0.435 
agrees with the H— and ZT-band ratios (Sluse et al. 



2008b 



|Fadely & Keeton|2011| >, although the latter two ratios differ by 
0.16 mag, suggesting that the //-band flux is still slightly mi- 
crolensed. At larger wavelengths, Fadely & Keeton (2011 1 find 
a //-band flux ratio B/A <0.164, suggesting that the H— and K— 
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band continua are microlensed as well, but probably by a mas- 
sive substructure. 

Broad Lines: They are apparently unaffected by microlensing 
unless M is significantly underestimated as suggested by the lit- 
erature data. If this scenario is correct, the broad Mg II and C III] 
emission lines should then be significantly microlensed. 

(f) FBQ 0951+2635: Deconvolution: The deconvolution is 
good. The residual flux in the vicinity of the lensed images 
amounts less than 0.03% of the brightest lensed image A and less 
than 0.2% of the flux of B. The lens galaxy is only ~0.15" ~ 4 
pixels away from image B and is only ~ 5 times fainter in the 
reddest part of the spectrum (i.e. > 7000A). This might lead to 
residual contamination of B by the lensing galaxy. The spectral 
regions 7200-7400 A and 7700-7950 A are unreliable due to an 
enhanced level of noise. 

Chromaticity: Despite the flat continuum flux ratio, we may not 
exclude possible residual contamination of image B by the lens 
because of the small separation between G and B. We identify 
ML of the continuum based only on the Mg II emission, which 
is not microlensed. The decomposition of the spectra around 
that line is however not entirely satisfactory because Fmh un- 
der Mgll does not show the monotonic variation expected if 
only the power law continuum was microlensed. This might be 
associated to microlensing of Fe II or to spurious effect of con- 
tamination by the lens. 

Notes: Microlensing of the continuum in this system, at differ- 
ent epochs, is supported by several other studies (e.g.|Schechter| 
let al. II 1998 1 |Jakobsson et al.||2005l |Munoz et al.||201 11 and ref- 



erence therein). A low amplitude chromatic change of B/A from 
4000 to 9000 A has been detect ed based on HST images ob- 
tained 2.5 years before our data (Mun oz et al.|201 1) . There is a 
good agreement between our spectral-based estimate of M and 
the radio flux ratio (which unfortunately lacks error estimates). 
The //-band ratio (Falco et a l.|1999[ l is larger by about 0.2 mag 
compared to our estimate. This offset is hardly explained by in- 
trinsic variability (Tab. C. 1 1 and therefore suggests that the con- 



tinuum is still microlensed in //-band. 

(g) BRI 0952-0115: Deconvolution: The deconvolution is 
very good. Spatially resolved narrow Ly a is visible in the back- 
ground image. This emission is not produced in the QSO but in 
the host galaxy of this remote quasar. 

Chromaticity: Because of the Lyman break, we do not estimate A 
in both sides of the emission but only in the continuum redward 
of Ly a. The ratio B/A seems however flat from the blue to the 
red with only an imprint of the broad emission lines. Therefore, 
the continuum emission is microlensed. 

Broad Lines: Our decomposition unveils ML of a significant 
fraction of the Ly a line. Because the flux leading to the absorp- 
tion is microlensed as the continuum, the latter does not appear 
in Fm, which unveils a (nearly) symmetric emission roughly 
centered on the narrow emission. 

Notes: The //-band flux ratio is similar to our continuum flux 
ratio, and therefore supports a significant microlensing (~0.45 
mag) at that wavelength. This is not a surprise as //-band corre- 
sponds to rest-frame UV emission {A ~ 3300 A,) which is small 
enough to be significantly microlensed. 

(h) SDSS Jl 138+0314: Deconvolution: The deconvolution 
is good with some residual flux and background excess in the 
vicinity of image B. Its origin is possibly associated to the QSO 
host galaxy. 

Chromaticity: There is a small chromatic change of C/B from the 
blue to the red part of the spectrum which is hardly explained by 
contamination from the lensing galaxy. The amplitude of this ef- 



fect is however very small and our measurement are compatible 
with no CML and no DE. 

Broad Lines: The ML of the BLR is large in this system, iso- 
lating the narrow component of the CIV flux in the non mi- 
crolensed fraction of the spectrum. The signal is less pronounced 
in C III]. We do not detect ML of the [He II] and [OIII] A 1663 
emission. If we use the Zf-band flux ratio as the correct value of 
M, we find significant microlensing of the C III] line. 
Notes: Our spectroscopic estimate of C/B differs by 0.25 mag 
from the Zf-band measurement (Blackburne et al. 2011). This 
is hardly explained by intrinsic variability and suggests signifi- 
cant reddening of image C. This also explains the //-band ra- 
tio. Alternatively, we might have underestimated the amount of 
microlensing in our spectra. This has to be confirmed with addi- 
tional data. 

(i) J1226-0006: Deconvolution: During the deconvolution 
process, we forced the separation between the 2 lensed im- 
ages and the lensing galaxy to match the HST separation in or- 
der to reduce cross-contamination. Although a good deconvo- 
lution is obtained, we keep in mind that residual contamination 
shouldn't be excluded due to the small separation of the system 
(A AG=0.437", A AB= 1.376") and of a pixel size twice larger 
than for the other lenses (i.e. 0.2"/pix.). 

Chromaticity: There is a strong change of the flux ratio from 
the blue to the red. We find evidence for unproper deblending 
of the QSO and of the lens flux, especially below A < 6500A 
-corresponding to the 4000A break of the lens- where the con- 
tinuum spectrum of the lens and of the QSO have similar shapes. 
The decomposition of the Mg II line also shows imprint of the 
H and K absorption bands from the lens. This effect could 
modify the intrinsic shape of the microlensing signal. Based 
on the Mgll and [Oil] ,13727 A lines, we find that only DB^J 
(M(blue) =0.850, M(red) = 0.765, fx =0.536) may cause the ob- 
served chromatic change, with image B being more reddened by 
the lensing galaxy. This is unexpected as the lens is closer from 
image A than B. Another possibility would be a color change 
associated to the intrinsic variability, as the time delay in system 
should be of the order of 25 days. 

Broad Lines: We observe ML of the blue component of Mg II but 
we are unsure of the role of the contamination by the lensing 
galaxy in our spectral decomposition. 

Notes: The //-band flux ratio is similar to the optical flux ra- 
tio A = 0.456 but deviates significantly from our line-based es- 
timates of M — 0.80. Although intrinsic variability and differ- 
ential extinction may play a role in the explaining the discrep- 
ancy, it seems plausible that the //-band continuum is in fact 
microlensed nearly at the same level as the optical one. 

(j) SDSS J1335+0118: Deconvolution: During the decon- 
volution process, we forced the separation between the bright 
lens image A and the lens galaxy component G to be identical to 
the HST separation. Small residual flux left after deconvolution 
close to A and B amounts less than 0.03% of the QSO flux. 
Chromaticity: The change of B/A from the blue to the red 
is incompatible with CML, as the latter has to be stronger at 
bluer wavelengths. Instead, we hypothesise DE (M(blue) =0.21, 
M(red) =0.23, fi =1.396), image A being more reddened than 
image B. We discard the possibility that the observed chro- 
matic trend is caused by residual contamination from the lensing 
galaxy as we estimate that more than 70% of the observed flux 
of the galaxy should contaminate image B to flatten the spectral 



10 Assuming the same M over the whole wavelength range implies 
that ML is the lower in the blue, in disagreement with the smaller source 
size at that wavelength. 
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ratio. 



Broad Lines: We observe ML of the red wing of C III] and of 
Mg II. For Mg II and C III] , there is a second (very)-broad com- 
ponent which is microlensed. 



Notes: The photometry published by Oguri et al. ( 2004 1, associ 



ated to data obtained about 2 years before our spectra, shows a 
slow decrease by 0.3 mag of Aiuba from g— to Zf -band, compat- 
ible with our results, but they argue this is not conclusive due to 
their photometric error bars. On the other hand, they find a flat 
spectral ratio B/A without clear imprint of the emission lines. 
This contrasts with our higher signal to noise spectra were we 
observe differential microlensing between the continuum and the 
emission lines. Owing to the expected IV (Table |CT| >, our esti- 
mate of M = 0.23 is compatible with the M = 0.29 measured 
in //-band. It however deviates significantly from the Zf -band 
measurement M = 0.41 (Ogu ri et al.|2004) >. To explain these ra- 
tios we have to postulate significant differential extinction. Oguri 
et al. (2004) observed a chromatic decrease of Am&i by 0.15 mag 
from r-band to Zf-band. This amount would lead to a value of 
M corrected from reddening M - 0.264. This still disagree with 
the Zf-band value but the two values get marginally compati- 
ble provided the effect of intrinsic variability is 50% larger than 
predictions from the structure function. 

(k) Q1355-225^) Deconvolution: During the deconvolu- 
tion process, we forced the separation between the bright lens 
image A and the lens galaxy component G to be identical to the 
HST separation. The deconvolution is good but low level resid- 
ual flux, up to 0.2% of the faintest image, is visible aside the 
faintest lensed image and the lensing galaxy. A small excess of 
flux appears in the background and in the PSF component of the 
galaxy at the wavelength of the peak of the Mg II emission. This 
flux is likely associated to image B but amounts less than 0.5% 
of the Mg II flux in that image. 

Chromaticity: Contrary to what is observed for the other sys- 
tems, the spectral ratio B/A in the continuum does not vary in 
a monotonic way (Fig. [CTTJc). The factor A is roughly the same 
in the continuum for A > 5450A but decreases significantly for 
bluer wavelengths. This trend cannot be explained by contami- 
nation from the lensing galaxy. There is significant differences 
between the spectra of images A & B at basically every wave- 
length suggesting a complex ML of the continuum and of the 
broad line region, including the region emitting Fell. Because 
we derive a similar value of M around Mg II and C III] , it seems 
plausible that the chromatic change of the flux ratio in the con- 
tinuum is caused by ML rather than DE although we cannot rule 
out the influence of the latter. The measurement of M from the 
[Ne V] narrow emission lines, although more uncertain, is com- 
patible with the one obtained for Mg II. 

Broad Lines: We observe ML of the red wing of Mgll. The 
C III] line being close to the red-edge of the spectrum, we es- 
timate A for this line in the range 4650-4680 A. We do not find 
ML of that line. 

Notes: Our estimate of M is in rough agreement with the 
/if -band flux ratio once we account of the possible effect of 
intrinsic variability (Table C.l i. The difference of ~ 0.15 mag 



between the H- and Zf-band ratios seems to be too large to 
be caused by intrinsic variability and is compatible with the 
//-band continuum affected by a small amount of microlens- 
ing. We scanned the spectra published by Morgan et al. (2003) 
and applied our decomposition method to these ones. We find 



forMgll (M,fi) = (0.33,0.63) and (M,/i) = (0.38,0.56) around 
CIII]. Like in our spectra, only the red wing of Mgll is mi- 
crolensed but not CIII]. The small differences on the derived 
values of M are easily explained by intrinsic variability. 

(1) WFI 2033-4723: Deconvolution: The deconvolution is 
good. Residual flux under the point-like images is < 0.1% of 
the QSO flux. 

Chromaticity: The amplitude of the chromatic differences be- 
tween the continuum in image B and C is small and probably 
caused by small uncertainties in the deblending of the QSO im- 
ages and of the lens galaxy. 

Broad Lines: There is ML of the blue wing of C III] or/and of 
Si III] but no ML of the Mg II line. 

Notes: Our estimate of M is in excellent agreement with 
the //-band and Zf-band measurements (Vuis soz et al.| [2008; 
|Blackburne et al.|201 1) >. 

(m) HE 2149-2745: Deconvolution: The deconvolution is 
good. Residual flux under the point-like images is < 0.09% of 
the QSO flux. The background flux retrieved by the deconvolu- 
tion process is probably associated to unproperly subtracted sky. 
Chromaticity: The measurement of A at the level of the C IV line 
is difficult because of its vicinity from the blue edge of the spec- 
trum. Our estimate of A for this line is performed in the range 
5070-5 120 A where Fell emission is minimal. There is no dif- 
ferential ML between the continuum and the emission lines. 
However the two spectra do not superimpose once scaled with 
the same magnification factor. First, there is a chromatic effect 
which may be caused by a small amount of DE, image B being 
more extinguished than image A, in agreement with its location 
closer to the lens galaxy. Intrinsic variability combined with the 
time delay of ~ 103 days (B urud et al.| 2002) might also explain 
the observed color difference. Second, the absorbed fraction of 
the CIV emission do not superimpose once scaled by B/A ~ 
0.242. This is likely caused by time-variable broad absorption 
which is seen in images A and B at two different epochs sepa- 
rated by the time delay. 



Notes: Comparison of our spectra with those of Burud et al. 
(2002) confirm our estimate of M — 0.245 at the level of the 
C III] line. This value is also in agreement with the //-band flux 
ratio. 



Appendix D: Characteristics of the extended 
sample 

We summarize here the main characteristics of the extended 



sample of objects introduced in Sect. 2.3 We provide the spec- 
tral ratio between the image pairs in Fig. D.l For each object, 
we also discuss our estimate of M from the spectra at the light of 
literature data (except for HE 2149-2745 which was discussed in 
Appendix O. 

(n) SDSS J0924+0219 (zi = 0.394, z s = 1.524): The VLT- 
FORS spectra of images A and B, obtained in Jan. and Feb. 2005 



by |Eigenbrod et al. ( 2006a i, show only a weak differential ML 



11 We found an error in the wavelength calibration of image B and 
therefore, we had to shift the original wavelength solution by 2 A for 
that image. 



between the continuum and the emission lines (C III] and Mg II) 
but no clear ML deformation of the lines (Fig. [4}. The spectral 
ratio is flat over the spectral range of our spectra. On the other 
hand, Keeton et al. (2006) identified flux ratios in disagreement 
with macro model predictions in both the continuum and emis- 
sion lines (Ly a, C IV and C III]) based on low resolution spec- 
tra obtained with the Advanced Camera for Survey (ACS) on 
29.05.2005. They demonstrated, based on microlensing simula- 
tions, that the observed anomalies in the image pair A- D could 
easily be explained by ML. Unfortunately, their spectra have in- 
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Object 


z 


At (days) 


Am (SF) 


(a) HE 0047-1756 


1.678 


(12) 


0.073 


(b) Q0142-100 


2.719 
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0.354 


(c) i>U!Sb JUZ4o-(J8z:> 


1.689 
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0.053 


(d) HE 0435-1223 


1.693 


6.5 


0.058 


(e) ^Ddb JUoUo+zUUo 
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(50) 


0.188 


(n) bUc>^ JtjVz4+U2iy 
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(5.7) 


0.061 


/ r-\ T7T3 /~\ f\c\ c i i i/;tc 

(1) rBQ 0951+2635 
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0.072 


(g) BK1 0952-01 15 


4AZ0 
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nil/; 
0.1 16 


(o) Jl 131-1231 


0.657 


12 T 


0.07 


(h) SDSS Jl 138+0314 


2.438 


(5) 


0.067 


(\\ n lift nnnA 


1 19^ 
1 . 1ZJ 


(10) 


VJ.lZj 


(j)SDSS J1335+0118 


1.570 


(49) 


0.162 


(k) Q1355-2257 


1.370 


(73) 


0.194 


(p)H1413+117 


2.55 


(20) 


0.116 


(1) WFI 2033-4723 


1.662 


62 


0.215 


(m) HE 2149-2745 


2.033 


103 


0.273 



Notes: t: Preliminary time delay estimate between A & B by Morg arTet al.| 
(2006) . 

Table C.l. Time delays for the observed systems and typical 
flux variation expected during this time-spent. When time de- 
lays have not been observed, we give in parentheses the pre- 
dicted value from our lens model. We do not quote the values 
for Q2237+0305 and for the image pair B-C in Jl 131-1231 be- 
cause the time delay is likely <1 day. 



sufficient spectral resolution to allow a proper investigation of 
the line profile differences. 

Notes: The fl ux ratio A/B = 0.44 + 0.04 in i/-band ( |Eigenbrod| 
et al. 2006a| l is in good agreement with our estimate of M. 



Although ^-band observation of this system exist ( |Faure et aL 
|201 l| l, they were not able to derive accurate photometry of the 
lensed images due to the ring. No other NIR/MIR photometric 
data of this system are available in the literature. 

(o) J1131-1231 (z, = 0.295, z s = 0.657): |Sluse^eTaTJ ( |2007l l 
presented single epoch spectra of this system obtained in April 
2003. They identified ML of the broad emission lines for im- 
ages A and C and presented MmD of the H B and Mg II emission 
lines using the image pairs A - B and C - B. They also identi- 
fied contamination of the spectra by flux from the host galaxy 
and empirically corrected for this effect. In image C, only the 
core of the emission line is not microlensed while microlensing 
affects the broad component of the lines in image A. Because of 
the uncertainty on the host contamination, exact values of M and 
fi may be more prone to systematic uncertaintie^than in other 
systems. Nevertheless, the microlensing of the emission lines is 
a robust result which qualitatively remains even if M or [i are 
under/over-estimated by up to 40%. 

Notes: The flux ratios in /T-band disagree with those derived 
from the narrow [O III] emission lines ( Slus e et al.|2007| [Sugai 
et al.|2 007 ). Differential extinction is not a plausible explanation 



because of the lack of monotonic change of the flux ratios with 
wavelength and because of the lack of strong hydrogen absorp- 
tion in X-ray (Cha rtas et al.|20 09). Therefore, it is likely that the 
TT-band flux ratio is significantly microlensed. 

(p) H1413+117 (zi unknown, z s = 2.55): This system is the 
first broad absorption line (BAL) quasars where ML has been 



unambiguously observed ( |Angonin et al. 1990 Hutsemekers 
1993). |Hutsemekers et al.| ( |2010| i presented the MmD for this 



system at four different epochs spanning a 16-years time-range. 
They identified that ML was affecting mostly image D, roughly 
in the same way, along the time-spent of the observations. 
The Lyo-, CIV, HB, and Ha emission have been analysed. In 
addition to the ML in the broad absorption, they found evidence 
for ML of the central core of the CIV and of the Lya (i.e. 
the wings are not microlensed), but no ML in the Balmer lines 
(see their figures 6, 7, 8). We show in Fig.|4]the decomposition 
for the CIV and for the HB lines for the spectra obtained in 
2005. Chromatic changes of the flux ratios are observed in this 
system ( |Hutsemekers et aL]2010| [Mufioz et al.|201l) . On one 
hand, there is differential extinction between images A & B, 
and on the other hand, there is chromatic microlensing of image 
D consistent with microlensing of a standard accretion disk 
( Hutsemekers et aT7||2.Q10|>. Note tha t microlensing also affects 
component C ( |Popovic & Chartas 2005), but the line profile 
differences are more subtle (jHutsemeke rset al.|2010| ). 

(q) HE 2149-2745 (z, = 0.603, z s = 2.033): Our spectra 
of this system do not show evidence of ML, however, Burud 
et al.| ( |2002| l presented spectra of the two lensed images of 
this BAL quasar, obtained on 19.11.2000, where they observed 
chromatic changes of the continuum. They mentioned subtle 
differences in the line profile of C III] but we fail to detect clear 
ML signature of the lines using the MmD on these spectra. The 
MmD is however difficult to perform because of the significant 
Fe II emission blueward of C III] and of the significant chro- 
maticity of the spectral ratio which is sensible even on the small 
wavelength range covered by the line. We display in Fig. [4] the 
decomposition of CIII] using fi(A) instead of the average ju 
between the blue and red part of the line. The chromatic changes 
observed in these spectra have been discussed by Burud et al. 
(2002) as possibly caused by differential extinction or chromatic 
microlensing. Although we may not rule out that differential 
extinction is present, we are now able to say based on the new 



spectra that the slope in the spectra of Burud et al. ( |2002| ) was 
mostly caused by CML. 

(r) Q2237+0305 (z, = 0.0394, z s = 1.695): There is a clear 
ML of the broadest component of the CIV and CIII] lines 
which has been observed in image A over the 3 years time-spent 



of the spectrophotometric monitoring presented in Eigenbrod 



et al. ( 2008 ). This signal has been used by Sluse et al. ( 201 1 1 to 
derive a size of the BLR in agreement with the size-luminosity 
relation obtained by reverberation mapping studies for other sys- 
tems. Image D has been found to be affected by differential ex- 
tinction. The MmD applied to the extinction corrected spectra, 
averaged over the first year of the monitoring (Oct. 2004-Sept. 
2005), is shown in Fig. [4] 

Notes: Although we originally used the lens model flux ratio 
A/D = 1.0 to make the MmD, we are also able to derive M em- 
pirically using the MmD. Following that procedure, we derive a 
very similar value of M (TableQ, but we disfavour M < 0.87 be- 
cause they lead to the appearance of a clear dip, that we consider 
as unphysical, in the center of the F Mfl fraction of the C III] line. 
Falco et al | (fT996) publis hed radio flux ratios A/D =0.77+0.23 
and |Minezaki et al.| ( p009l ) published D/A = 0.87 + 0.05 at 1 1.67 
fim. These values are in good agreement with our spectroscopic 
estimates. We should however notice that small differences be- 
tween the MIR flux rat i os of |Minezaki et aL] ( |2009[ ), |Agol et al.] 
(2009) and Agol et al. (2001|), as well as chromatic changes in 



12 Because the narrow [OIII] lines are resolved and possibly macro- 
magnified by a different amount than the broad lines, their use to set-up 
M might be another source of systematic errors. 



the MIR, suggest that a small amount of microlensing may still 
affect these wavelengths. 
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